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Abstract 

We investigate the cosmological dynamics of a brane Universe when 
quantum corrections from vacuum polarization are taken into account. 

New vacuum de Sitter points existing on Randall-Sundrum brane are de¬ 
scribed. We show also that quantum correction can destroy the DGP de 
Sitter solution on induced gravity brane. 
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Investigations of quantum effects in a strong gravitational field and their 
applications in cosmology have a long story. Since the beginning of 70-th many 
interesting results modifying the standard Friedmann cosmology due to vac¬ 
uum polarization and particle production have been obtained. In particular, a 
vacuum de Sitter solution and an inflationary regime driven solely by vacuum 
polarization without any matter jTj was described even earlier than a common 
scalar field inflationary scenario. A detailed analysis of cosmology with vacuum 
polarization (we will consider only this effect in the present paper) have been 
done in (2. After some period of stagnation, this problem begins to attract a 
considerable attention last several years, mainly due to development of modified 
gravity models. The form of extra terms in cosmological equations of motion 
caused by vacuum polarization does not depend on a particular theory of grav¬ 
ity, however, peculiarities of a background metric in these theories could result 
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in some new dynamical regimes. Recently modifications caused by the vacuum 
polarization have been studied for regimes with soft future singularities. Such 
regimes, being impossible in the standard cosmology, are rather typical in some 
modern cosmological scenarios (see for example, 0), in particular, they are 
present in induced gravity brane models EJ. Quantum corrections change the 
dynamics significantly, leading to a softer singularity or even to non-singular 
solutions Emm. 

Another interesting problem is stability of classical solutions with respect 
to quantum corrections. As quantum terms contain higher time derivatives 
in comparison with corresponding classical equation of motion, some classical 
solutions may become unstable. In jTj this instability is described for certain 
regimes in induced gravity brane cosmology. 

It is well known that the vacuum polarization leads to the following vacuum 
expectation value for the energy density: 

p q =< T 0 o >= k 2 H A + k 3 {2HH + 6 HH 2 - H 2 ), (1) 

where H is the Hubble parameter, k 2 ,k 3 depend upon the spin weight of the 
different fields contributing to the vacuum polarization. 

If the classical Friedmann equation has the form of algebraic dependence of 
the Hubble parameter upon the energy density of the Universe H = H(p) 1 sub¬ 
stituting p —> p+p q we get a differential equation which governs the cosmological 
evolution with quantum corrections. 

It is clear that k 2 - term in p q can be incorporated into the function H(p). 
Only ^ 3 -term containing time derivatives may provide instability. In the paper 
(7 the case of k 2 = 0 have been studied in braneworld models. It is a good 
assumption near a soft future singularity where H —> oo while H is finite, and 
fc 3 -term in (1) dominates. In the present paper we describe general features of 
brane dynamics with quantum corrections in their full form. 


2 


We start with modification of classical equations due to fo-term. After 
that we consider the problem of stability. Plots of the function H(p) (which 
includes fe-ff 4 term) presented below can be understood in two different ways. 
First, the curve H(p) can be considered as a track of cosmological evolution in 
some specific regime when ks-teim can be neglected. In this interpretation, p 
is the sum of matter density p m and brane tension a. The evolution runs from 
higher to lower energy (from the right to the left in the plots) till the point 
Pm = 0, p = a is reached. On the other hand, these plots can be interpreted 
as sets of de Sitter fixed points (with p = a and H = 0) of a general dynamics 
with non-zero k%. 

Before we discuss branes it is useful to remember known results in the stan¬ 
dard scenario. The modification of the standard Friedmann cosmology caused 
by the ^ 2 -term is shown in Fig.l. The classical cosmology corresponds to the 
straight line H 2 ~ p. The case < 0 results only in changing this dependence 
to H 2 ~ y/p for large H. A positive &2 modifies the situation significantly. In 
particular, we can see a new vacuum de Sitter point [ 8 ]. This point is unstable if 
&3 < 0, on the other hand, the point (0,0), representing the late-time Friedmann 
regime is stable for /C 3 < 0. It means that a trajectory starting in the vicinity of 
the vacuum de Sitter point has inflationary behavior at the initial stage, then 
it leaves a neighborhood of de Sitter due to instability. As a result, inflation 
ends, and finally the trajectory reaches the Friedmann late-time attractor. This 
scenario of acceleration expansion with the natural exit without any kind of 
classical matter is called as ’’Starobinsky inflation” [I]. 

We now turn to braneworld models. In order to understand how k %-term 
modifies the Randall-Sundrum (RS) brane it is necessary to remember about 
two branches, arising from p 2 term in the equation of motion. The unmodified 
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Figure 1: The p{H) dependence in the standard cosmology without &2 term of 
the quantum corrections (solid), with negative Ji 2 (dashed) and with positive Jt 2 
(dot-dashed) 

equation for the RS brane is 

H 2 = A/6 + p 2 /{9M 3 ). (2) 

Here A is a cosmological constant in the bulk, which should be negative in a 
realistic model, M is the 5-dimensional Planck mass. The equation (2) when 
solved with respect to p has two solutions for a given H, one with (+)-sign and 
the second with (—)-sign before the square root. This equation is represented as 
a shifted parabola on the plane (H 2 , p). As the negative brane tension ultimately 
leads to instability, the second branch of the parabola (which corresponds to 
negative p) is unphysical. Positive k 2 -term changes this classical picture in a 
way, very similar to the standard cosmology case (see the dashed curve in Fig.2) 
with the (—) branch remaining the unphysical one. A vacuum de Sitter point 
now corresponds to a positive root of the forth-order equation 

(k 2 /9M 3 )H 8 — H 2 + A/6 = 0. (3) 

In the limit At 2 —► 0 the value of Hubble parameter in this point tends to infinity. 
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Figure 2: The Randall-Sundrum cosmology without /c 2 -term (solid), with neg¬ 
ative &2 (dot-dashed) and with positive &2 (dashed) 

The existence of a vacuum solution in the RS brane theory with the quantum 
corrections was first noticed by Nojiri and Odintsov in Emni 

The picture for negative &2 is shown in the dot-dashed curve in Fig.2. The 
(+) branch is similar to the standard cosmology, however, the (—) branch is also 
shifted partly to a physical domain p > 0. In particular, it has a new vacuum 
de Sitter point. The location the this point is also given by a positive root of 
(3), however we will see below that stability properties of these two de Sitter 
points are different. 

Finally, we consider induced gravity (IG) brane. The classical equation of 
motion is 

m 4 (ff 2 - p/(3m 2 )) 2 = M 6 (H 2 - A/6). (4) 

Here m is the 4-dimensional Planck mass. As the equation of motion is quadratic 
with respect to p there are also two branches of solutions. It should be however 
noted that in some situations a negative brane tension on an IG brane does 
not ultimately lead to instability EH, so negative p may have a physical sense 
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in contrast to a RS brane. The plots in Fig.3 have been done for the case of 
Minkowski bulk (A = 0). The (+) branch of the solid curve has a Friedmann- 
like limit when the energy density of the brane vanishes ( H —> 0 if p —> 0), the 
(—) brane has a vacuum de Sitter point, founded first by Dvali, Gabadadze and 
Porrati in [21 ■ We will call it DGP point. 

A negative (the dot-dashed curve in Fig.3) does not change the config¬ 
urations of branches. On the other hand, positive &2 can modify the picture 
seriously. Large enough positive &2 transforms the diagram so that the lower 
branch never enters into the p > 0 half-plane, and the general picture resembles 
the case of RS brane (compare the short-dashed curve in Fig.3 with dashed 
curve in Fig.2). It is evident that the DGP point is absent in this case. In 
can be easily calculated that this situation is realized if &2 > (4/9)(to 6 /M 6 ). 
Nonzero A in the bulk does not change this general picture qualitatively, though 
the expression for the critical &2 becomes less simple. We can see also that the 
upper branch has a vacuum de Sitter point for any positive k 2 - 

As it was shown recently in (7j , the effective phantom branch (when dH/dp 
is negative) of IG brane is unstable with respect to fc 3 -term in the vacuum 
polarization. Here we will show that this result is rather general and does not 
depend on a particular kind of a brane as well as on the value of k- 2 - 

We consider a general situation when matter density in the corresponding 
effective Friedmann equation is some algebraic (may be multi-value) function of 
the Hubble parameter H , so the evolution equation for a brane has the form: 

p + a = F i (H), (5) 

where the index i marks different branches, a is the brane tension. It is clear 
that both RS and IG branes are particular cases of (5). 

In the de Sitter point the matter density is diluted, so the equation becomes 
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Figure 3: The induced gravity brane cosmology without A/ 2 -terrn (solid) , with 
negative k? (dot-dashed), with positive &2 which allows the DGP point, denoted 
as A (long-dashed), and with positive &2 which does not allow the DGP point 
(short-dashed) 


a = Fi(H). 

Now we add the quantum corrections (1). The fe-term can be incorporated 
into the algebraic functions 

F i {H) = F i (H) + k 2 H 4 . 

The new de Sitter point is given by 

a = Ft(H). 

The resulting second order differential equation for if can be written in the 
form of a system of two first order equations 
H = C, 

C = -3 CH + g + 2WI^(Fi -a) = MH, C ). 

There is an equilibrium point of this system (if, C ) = (0,0), which corresponds 
to the de Sitter solution. In order to investigate stability of this equilibrium 
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point we need a linearized system: 


C=(%)oC + (§^) 0 H, 

H = C. 


The eigenvalues of this linearized system have the form: 


Ml ’ 2 = 2 


y dC’ 


( 9f_\ +J(2L\2 + a(^L\ 

V a/-<) 0 ^ \/\ an )0 ^lauJo 


dC' 


dH' 


We can easily see that the first eigenvalue is negative if and only if ( : §jj)o < 0. 
The second eigenvalue is always negative (because (§q)o = —3 H is negative in 
an expanding universe). 

Let us now evaluate 


A - 1 

y dH >0 2H 2 k 3 


—<j + FA 


1 dF 


2Hk 3 dH 

In the DeSitter stable point —a + Fi = 0, and we have finally 


( ?A) 1 

[ 8H >0 2 Hk 3 dH l 

As a result, if > 0 (a normal branch), the de Sitter solution is stable if 
k 3 < 0 (as in the standard cosmology). In the opposite case (a phantom branch) 
we need k 3 > 0 for stability. 

We can conclude, that if we have k 3 < 0 in our Universe (which is needed 
for stability of the Minkowski space), all de Sitter points on phantom branches 
are unstable. 

Note, that we did not specify any particular modified gravity theory, which 
may even be inspired by some scenario different from brane models. All we need 
is a function H(p). Of course, this analysis can not be applied to a situation 
when equation which relates H and p is a differential one (the most important 
example of this case is a scalar field with non-zero potential playing the role of 
matter). 

It also should be noted, that, strictly spiking, we have proved the following 
statement: A de Sitter point, being a future attractor of a Universe is unstable 



with respect to vacuum polarization if it located on a phantom branch. However, 
the results of [3 indicate that an effective phantom regime is unreachable during 
a cosmological evolution of a IG brane when quantum corrections are taken into 
account, and it is quite reasonable to suggest that quantum corrections prevent 
a classical phantom regime from realization in a general situation. This problem 
requires further investigations. 

In any cases, we have obtained new restrictions for possible explanations of 
a present phantom-like state of ” dark energy” which is often claimed as being 
favorable by observations ismansi. Several years ago it was remarked that a 
matter with a constant equation of state parameter u> < — 1 causes an unwanted 
Big Rip future singularity HM3IH]. This singularity can be avoided in some 
scenarios with a phantom scalar field (a scalar field with the wrong sign of the 
kinetic term) O 2} , however, the fact that the energy of a phantom scalar 
field is unbounded below causes severe problem in quantum theory EH ( recently 
proposed more complicated models for a scalar field phantom see, for example, 
in mmi). All these problems are absent in modified gravity proposal. In 
this approach the matter in the Universe remains standard (and, so, there are 
no problems with matter instabilities), and the effective phantom behavior is 
achieved due to significant modification of the Friedmann equation. The IG 
brane is a famous example of such kind of theory E2 Our results, showing 
that quantum corrections have significant influence exactly on those branches 
of cosmological equations which simulates a phantom behavior, make clear that 
the modify gravity proposal is also not free from instability problems. 

Remembering shapes of H(p ) dependence in braneworld models, we can see 
that the new branch of RS brane existing for k -2 < 0 is stable, as well as new 
vacuum de Sitter point (3) on this branch. This point has no analog in standard 
cosmology and is similar to DGP point on IG brane. On the other hand, this 
point in the case of &2 > 0 is unstable (like in standard cosmology), giving a 
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possibility of realization of Starobinsky inflation on RS brane. We also can easily 
see that the Starobinsky inflation on IG brane is always possible (of course, if 
&2 > 0) on the (+) branch and for &2 < (4/9)(m 6 /Af 6 ) on the (—) branch. 

We have studied modifications of brane cosmology caused by vacuum polar¬ 
ization. This effect consists of two different terms in an effective energy density. 
The term proportional to H 4 alters possible fixed points of a cosmological dy¬ 
namics, while the term containing time derivatives of H may change stability 
properties of these fixed points. A general condition for the future fixed de Sit¬ 
ter point to be stable have been derived. We should however note that all these 
results have been obtained when possible quantum corrections in the bulk are 
neglected. This assumption is reasonable in studies of quiescent future singular¬ 
ities on a brane, because they are singularities of embedding [3J, while the bulk 
remains regular (and, so, far from a quantum regime). Is quantum corrections 
in the bulk important for a brane dynamics in general situation remains unclear, 
we leave this problem for future investigations. 

Acknowledgments 

This work is supported by RFBR grant 05-02-17450 and scientific school grant 
2338.2003.2 of the Russian Ministry of Science and Technology. Authors are 
greatful to Alexey Starobinsky, Varun Sahni and Yuri Shtanov for discussions. 

References 

[1] A. Starobinsky, Phys. Lett. B91, 99 (1980) 

[2] M.Fischetti, J.Hartle and B. Hu, Phys. Rev. D20, 1757 (1979). 

[3] J. Barrow, Class. Quant. Grav., 21, L79 (2004) 

[4] Yu. Stanov and V. Sahni, Class. Quant. Grav., L101 (2002). 


10 


[5] S. Nojiri and S. Odintsov, Phys. Lett. B595, 1 (2004) 


[6] S. Nojiri, S. Odintsov and S. Tsujikawa, Phys. Rev. D71, 063004 (2005) 

[7] P.Tretyakov, A.Toporensky, Yu.Shtanov and V.Sahni, ’’Quantum effects, 
soft singularities and the fate of the universe in a braneworld cosmology”, 
gr-qc/0510104 

[8] S.Mamaev and V.Mostepanenko, Sov.Phys.JETP 51, 9 (1980). 

[9] S. Nojiri and S. Odintsov, Phys. Lett. B484, 119 (2000). 

[10] S. Nojiri and S. Odintsov, Grav. Cosmol. 8, 73 (2002). 

[11] Yu. Shtanov and A. Viznyuk, Mod.Phys.Lett. A20, 2283 (2005). 

[12] G. Dvali, G. Gabadadze and M. Porrati, Phys.Lett. B 485, 208 (2000) 

[13] U.Alam, V.Sahni, T.Saini and A.Starobinsky, Mon.Not.Roy. Astron.Soc., 
354, 275 (2004). 

[14] J. Alcaniz, Phys. Rev. D69, 083521 (2004). 

[15] T.Roy Choudhury and T.Padmanabhan, Astron. Astrophys. 429, 807 
(2005). 

[16] A. Starobinsky, Grav. Cosmol. 6, 157 (2000). 

[17] V. Sahni and A. Starobinsky, IJMPD 9, 373 (2000). 

[18] R, Caldwell, Phys. Lett. B545, 2329 (2002). 

[19] P.Singh, M.Sami and N.Dadhich, Phys.Rev.D68, 023522 (2003). 

[20] M.Sami and A.Toporensky, Mod.Phys.Lett. A19, 1509 (2004) 

[21] J.Cline, S.Jeon and G.D.Moore, Phys.Rev.D70, 043543 (2004) 


11 



[22] F.Carvalho and A.Saa, Phys. Rev. D70, 087302 (2004). 

[23] A.Andrianov, F.Cannata and A.Kamenshchik, ’’Phantom Universe from 
CPT symmetric QFT”, gr-qc/0512038 

[24] V. Sahni and Yu. Shtanov, JCAP 0311:014 (2003). 


12 



